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We show that it is possible to initialize and manipulate in a deterministic manner protected qubits 
using time varying Hamiltonians. Taking advantage of the symmetries of the system, we predict 
the effect of the noise during the initialization and manipulation. These predictions are in good 
agreement with numerical simulations. Our study shows that the topological protection remains 
efficient under realistic experimental conditions. 



Introduction. — Quantum error correcting codes are 
deemed a crucial tool for the implementation of quantum 
algorithms p] [2] ■ Soon after their discovery, it was shown 
that some errors can be self-corrected by assigning them 
an energetic penalty. This is based on the encoding of 
a single logical qubit in the degenerate ground state of 
a composite system whose excited states are separated 
by a finite energy gap from the ground state subspace 
[3H5] . Following the pioneering proposal by Kitaev [6] , it 
has been shown that when this scheme is implemented 
using specific Hamiltonians with collective symmetries, 
it can reduce exponentially the influence of decoherence 
on logical qubits and quantum gates [TMllj. 

It is usually recognized that logical qubits should verify 
the so-called di Vincenzo criteria [12]. In particular, it is 
necessary to be able to initialize to a fiducial state and 
perform single and two qubits rotations on them. This 
initialization is difficult because the protected states are 
highly entangled states which are thus not reachable from 
simple, separable states using local operations. Also, for 
a protected qubit, there is an apparent conflict between 
the protection and the ability to manipulate it. We show 
here how to overcome these difficulties using time vary- 
ing Hamiltonians with well-defined symmetries. Using 
both high-order perturbation theory and numerical sim- 
ulations, we also take into account explicitly the effect of 
random noise during all these processes showing that the 
qubits remain protected during these operations. 

Protection Hamiltonian. In Rcfs.[6HTT]. well- 

designed Hamiltonians acting on an array of N x N 
two-level systems were introduced and shown to yield 
degenerate quantum states largely immune to decoher- 
ence. The protection relies on the existence of a set 
of symmetries, which are collective operators acting on 
rows or columns independently, i.e. P\ = YljLi ®l p Qj = 
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IT—I of * where <rf, is the standard Pauli matrix in the 
a direction in spin-space acting on the two-level system 
located at the intersection of row i and column j. They 
obey the following rules : [Pi, Pi'] = [Qj,Qj'] = = 
{PiiQj}) Pi =1 = Qj- The symmetries enforce the two- 
fold degeneracy of all eigenstates, including the ground 
states {|0i), It is convenient to choose these states 

so that Pi acts as a logical Pauli operator denoted f z 
while Qj acts as f x . These states are complex entangled 
states composed of arrays of TV 2 two-level systems that 
cannot be prepared nor manipulated easily. They are 
separated from the other eigenstates by a finite energy 
gap A, whose exact value depends on the specific model 
and on the system size, lines More quantitatively, the de- 
generacy lifting when the system is subjected to a static 
noise of amplitude h independently acting on N of its 
constituents scales as h N /A N_1 , which decreases expo- 
nentially with N provided h is smaller than A. Accord- 
ingly, the coherence time increases exponentially with TV. 

In the physically relevant case of a dynamical noise, 
this exponential improvement holds only below a criti- 
cal size N* which can roughly be estimated as the ratio 
between the energy gap A and the typical energy scale 
{ksT for a noise source in thermal equilibrium) which 
the environment can provide to the system in a single 
elementary process. The existence of this finite N* is 
consistent with the notion that topological order in 2D 
systems disappears at finite temperature in the thermo- 
dynamical limit 13-15 . However, for the purpose of 
quantum computation, we are interested in small sys- 
tems initially prepared in their ground states and time 
scales much shorter than the thermalization times. 

To illustrate numerically how one can circumvent the 
difficulties of initializing and manipulating such pro- 
tected states, we will use the long- range Hamiltonian 

N I N \ 2 n / n \ 2 



2 



where J x and J y are numerical constants and which can 
be realized either using trapped ions [TT] or supercon- 
ducting circuits [Ttj] . 

Initialization. — We now turn to the state initializa- 
tion. In the same spirit as [T7], we proceed by apply- 
ing a static field that acts independently on each two- 
level system. Such a local field is then adiabatically 
turned off from t = to t = T while the protecting 
Hamiltonian is turned on during the same interval, in 
a process described by the time dependent Hamiltonian 
H t = f(t)H + (l-f(t))S v , where f(t) is a slowly varying 
function such that /(0) =0 and f(t) t ^t±£ i with a char- 
acteristic time r, and S v — J2i j a ij describes the effect 
of the static field. The system can be prepared through 
single site operations in the (ferromagnetic) ground state 
of S y . This state is an eigenstate of Pi with the eigen- 
value 1 which ensures that, after the evolution, the final 
state is the protected eigenstate with the same eigen- 
values, |0i). In such a state preparation scheme, the 
system only reaches asymptotically the decoherence pro- 
tected state. The time evolution of the instantaneous 
energy levels is indicated in Fig. [IJ (inset) : at t = the 
energy levels are those of S v and close to those of Hq 
at t = 5r (where we have taken f(t) = exp(— t 2 /r 2 )). 
The effect of this procedure on the logical state is shown 
in Fig. [IJ top where we have plotted the error 1 — F(t) 
where F(t) is the fidelity F(t) = \(0 L \ip(t))\ 2 of the pre- 
pared state, with \tp(t)) the quantum state of the system 
at time t. This figure shows that for r = 100, the final 
state is extremely close to the target state(error of the 
order of 10~ 10 ). Indeed, the error after a time t = 5r 
diminishes as t increases and the error is below 10~ 5 for 
r > 20 as shown in Fig. [TJ bottom. 

However, in a real experiment, during the whole evo- 
lution, the system is subjected to noise that will ren- 
der state preparation imperfect. An important point is 
to evaluate the effect of noise on the fidelity. Simula- 
tions show that the preparation scheme is quite robust 
to noise. We have plotted in Fig. [TJ the error for different 
noise amplitudes as a function of t. We consider here a 
static noise with a random orientation applied on each of 
the individual two-level systems. This figure shows that 
even for realistic error rates, on the order of 0.01J, there 
exists some values of t for which the error will remain 
smaller than 10~ 3 (2x2 case). 

Our results are important for the realization of experi- 
mental tests of protected qubits as they show that, under 
realistic experimental conditions, errors will not destroy 
the quality of the prepared state. Let us now turn our 
attention to the manipulation of a single protected qubit. 

Manipulation. — We address the possibility of manip- 
ulating the ground states of the protecting Hamiltonian 
using local operations. 
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FIG. 1: (Color online) Top: Example of the state evolution 
as a function of time without noise in 2x2 case and for r = 
100 which is large compared to A" 1 = (\/2 + l)/2. Inset : 
Evolution of the energy levels (normalized to J = J x = J y ) 
as a function of time: at t = the energy levels are those of 
S y and close to those of Ho at t = r. Bottom : Error as a 
function of time for a static noise applied on all qubits on a 
2x2. Noise amplitude 10~ 2 (squares) and 10 -3 (circles) 

We remark that, due to the symmetry of the system, an 
operator acts within the logical qubit space if it either 
preserves or changes into their opposites simultaneously 
all the eigenvalues of the symmetry operators. This con- 
dition will be denoted in the following Condition 1. 

Firstly, we study the noiseless case. We thus consider 
a perturbation Hamiltonian of the form S u = . afj 
where u = x or y. This Hamiltonian is assumed to be ap- 
plied using a smooth time- varying gate function so that 
the total Hamiltonian reads H t = H a + g u (t)S u . The 
maximal value attained by g u (t) , gm ax remains small 
compared to the energy gap of Hq so that no level cross- 
ing occur. Furthermore, the typical variation time is cho- 
sen much smaller than the inverse of the gap frequency 
so that the transformation can be considered as adia- 
batic within the excited manifold (it will naturally not 
be adiabatic within the protected state where the energy 
difference is initially zero). According to the quantum 
adiabatic theorem, these conditions are sufficient to en- 
sure that the final state will remain within the protected 
subspace (|0l),|1l)). Thus, the unitary operator U de- 
scribing the transformation from initial to final state can 
be decomposed over the logical Pauli operators, denoted 
r x ' y ' z , and the two-dimensional identity matrix, I 2 

u = ai i 2 + a ^ x ( 2 ) 

i=x,y,z 
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In order to study the effect of the manipulation, one can 
apply high order perturbation theory and decompose the 
result according to |2]). At the M th order, one obtains 
a sum of terms of the form H M = ^i u j^t 2 ,n ' ' ■ Cjm 
(u is an arbitrary direction). Due to the adiabatic theo- 
rem, Hm will act within the logical qubit space so that 
it must verify Condition (1) to have a non trivial effect. 
Two pairs of exclusive possibilities can then be consid- 
ered: 

1. the hrst pair concerning the Pi operators (whose eigen- 
values are denoted pi) 

(a) Hm conserves all pc in this case, Hm commutes 
with the symmetry operators Pi and thus with the logi- 
cal operator f 2 . It must then induce a term of the form 
aS-2 + bf z where, here and in the following, a and b are 
constants depending on the exact shape of g u ; 

(b) Hm changes all pi in their opposites: in this case, 
Hm anti-commutes with the symmetry operators Pi and 
thus with the logical operator f z . It must then induce a 
term of the form af x + bf v ; 

2. the second pair concerning the Qj operators (whose 
eigenvalues are denoted qj) 

(c) Hm conserves all qj: in this case, Hm commutes 
with the symmetry operators Qj and thus with the log- 
ical operator^. It must then induce a term of the form 
al 2 + bf x ; 

(d) Hm changes all qj in their opposites: in this case, 
Hm anti-commutes with the symmetry operators Qj and 
thus with the logical operator f x . It must then induce a 
term of the form ar y + bf z . 

For a given Hamiltonian Hm, only one term will re- 
main for the two pairs of conditions which can be sum- 
marized in Tab. We see that, for each of the four 
conditions, one term and only one will survive which sim- 
plifies greatly the analysis. 

Condition (c) Condition (d) 

(gj -» 9j > V? ) (ij -» -gj j Y?) 

Condition (a) » z 
(pi — > pi, Vi) 

Condition (b) „ x „ y 

(Pi ->■ -Pi.Vi) 

TABLE I: Effective operator induced by the manipulation as 
a function of its effect on the ground-state quantum numbers. 

As discussed above, high order perturbation theory 
yields a sum of terms of the form Hm,M € N. Each 
of these terms will have an effect summarized in Tab. |l| 
and we have to consider the dominant terms in the sum. 
Let us first consider the case of a manipulation along y. 
At any order of perturbation theory, S y will conserve the 
Pi quantum numbers, all terms appearing in the pertur- 
bation theory development are thus in the first line of 



Tab. |[]). For M < N, there is at least one qj quantum 
number which is preserved so that Condition (1) can not 
be fulfilled: the manipulation in this case has no effect, 
in good agreement with the idea of protection. The first, 
non-zero term will appear when all columns are touched 
at least once by the manipulation (condition (d)), that 
is M = N. In this case, the effect will be on the order 

of ( Sm £* c ) t z where t is the duration of the manip- 
ulation. A similar reasoning can be made for a manipu- 
lation along x which preserves the qj quantum numbers: 
again, the first, non-zero term will appear when all lines 
are touched at least once by the manipulation, that is 
M = N. In this case, the effect will be of the order 

of ( 9t 2^ x J ^t x . Finally, combining these two types of 

operations, any arbitrary rotation on the protected qubit 
can be performed, a key operation toward the use of such 
systems for quantum calculation purposes. 

Let us now tackle the more difficult question of the 
resistance of this procedure to noise. We consider here 
an arbitrary static noise which is expressed as a sum of 
single spin terms. In practice, one will overestimate the 
effect noise by taking the dominant noise assumed to be 
along v with amplitude The total Hamiltonian now 
reads H — H + g(t)S u + fS v . 

At the M th order of perturbation theory, one 
is left with a sum of terms of the form Hm = 

first term comes from the manipulation term and the sec- 
ond from the noise term. We have naturally m + n — M. 
In the case where manipulation and noise are parallel 
(u = x = v or u = y = v), the two terms add up to 
yield a term of the form (g^ ax + f) N r x ' z . In this case, 
the rotation axis is not modified but the rotation angle is 
affected to first order in f/g^ax^ without any gain from 
the protection. 

We now turn our attention to a case where manip- 
ulation and noise are not parallel, for instance u = y 
and v = x. The two terms affect respectively only the 
columns (manipulation) or the lines (noise). The previ- 
ous result can be used: one can have a change of sign 
of a family of quantum numbers only if the correspond- 
ing power (to for the manipulation and thus the columns 
and n for the noise and thus the lines) is larger than or 
equal to N: the smallest order for noise is thus either 
2 (when there is no sign change of the p^) or N (when 
all the pi change sign) while the smallest order for the 
manipulation is either (when there is no sign change of 
the qj) or N (when all the qj change sign). Finally, these 
terms must be compared with the effect of the manipula- 

tion (proportional to ( g "£ x ) ^) as it is the ratio which 
is relevant and the corresponding results are detailed in 
Tab. |ll]). Similar results can be obtained for the other 
cases where noise is not parallel to the manipulation. We 
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Condition (c) Condition (d) 

(qj ->• g 3 -,Vi) (qj-*-qj,Vj) 



Condition (a) 
(Pi -> Pi,Vi) 
Condition (b) 
(Pi -> -Pi,Vi) 



not relevant 



A2 ' 



TABLE II: Dominant term induced by the manipulation and 
noise terms as a function of its effect on the quantum numbers 
where amplitudes are normalized by the effect of the manip- 
ulation in the absence of noise. 




0.000 0.002 



0.004 0.006 
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FIG. 2: Deviations of \a x \ (red circles), lo^l (green squares) 
and \a z \ (blue rhombi) away from the noiseless manipulation 
case as a function of the error amplitude / for a 2 x 2 array. 
The manipulation amplitude is jJJ, al = 0.304 corresponding 
to a f phase-shift. 



have thus demonstrated that manipulation of the logical 
qubit state can be performed in a deterministic manner 
and that this procedure is relatively resistant to noise. In 
the worst case (noise parallel to the manipulation), the 
noise sensitivity is the same as in the absence of protec- 
tion, while for the other cases it varies at least quadrati- 
cally with the perturbation strength, and the deviations 
of the rotation axis in qubit space are exponentially small 
in N. 

In order to confirm these results, we used numerical 
simulations performed on 2 x 2 arrays. We first checked 
the effect of a manipulation along x in the absence of 
noise. As expected, the coefficients along a y and a z are 
close to zero (typically smaller than 10 -9 .), showing that 
the manipulation is a pure rotation along z of the log- 
ical qubit. ot\ and a x are, also with a good precision, 
sinusoidal functions of the manipulation amplitude, g x . 
This can be seen as Rabi-like oscillations between the 
two-qubit states. 

Concerning the effect of the noise, the effect of the pro- 
tection is most apparent for a noise perpendicular to the 
manipulation axis. Although we have explored numeri- 
cally all combinations of directions for the manipulation 
and noise, we discuss as above the specific case u = y 
and v = x (Fig. [2]). 

Without noise, the manipulation should then be equiv- 



alent to a perfect dephasing. We verify that the devia- 
tions of a x and a y away from a noiseless manipulation 
have a quadratic variation with the noise amplitude as 
expected since N — 2 (see Tab. dill), bottom row). Note 
that the deviation of a z also varies quadratically even 
though it is several orders of magnitude smaller. 

In any physical implementation, the binary interac- 
tions between spins in the protected Hamiltonians are 
likely to be subjected to residual fluctuations, so that 
they become space and time dependent. In contrast with 
early proposals of self-correcting qubits [3HS] , these fluc- 
tuations do not lift the ground state degeneracy and do 
not induce decoherence because they preserve the non- 
local symmetries P$ and Qj . This remains valid as long as 
the gap remains open, that is as long as the fluctuations 
in the coupling are small in relative value. During the 
manipulation, the effect of these fluctuations can be an- 
alyzed as discussed above. It appears immediately that 
they have no effect on the rotation axis. However, as in 
the case of a single spin noise parallel to the manipu- 
lation axis, the rotation angle is affected by an amount 
proportional to the noise amplitude. 

Conclusion. — Using high order perturbation theory 
and exploiting the large number of symmetries present in 
Hamiltonians yielding protected qubits, we have shown 
that these qubits can be initialized and manipulated in 
a noise resistant manner. This shows that there is no 
contradiction between the concept of protection and the 
ability to manipulate a qubit. 
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